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Abstract 

The pair production of Z bosons is studied using the data collected by the L3 
detector at LEP in 1998 in e"'"e~ collisions at a centre-of-mass energy of 189 GeV. 
All the visible final states are considered and the cross section of this process is 
measured to be 0.74lo;^4 (stat.) ib0.04 (syst.) pb. Final states containing b quarks 
are enhanced by a dedicated selection and their production cross section is found to 
be 0.18to;Q7 (stat.) ±0.02 (syst.) pb. Both results are in agreement with the Standard 
Model predictions. Limits on anomalous couplings between neutral gauge bosons 
are derived from these measurements. 
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1 Introduction 



Since 1997 LEP is running at centre-of-mass energies, ^/s, above the production threshold of 
Z boson pairs. This process is of particular interest as it constitutes an irreducible background 
for the search of the Standard Model Higgs boson and to several other processes predicted by 
theories beyond the Standard Model. In addition it allows the investigation of possible triple 
neutral gauge boson couplings, ZZZ and ZZ7 [T, forbidden at tree level in the Standard Model. 

The experimental investigation of ZZ production is made difficult by its rather low cross 
section, compared with competing processes that constitute large and sometimes irreducible 
backgrounds. The existence at threshold of this process was established at y/s = 183 GeV [2]. 
In the following, the analysis of the data collected at 189 GeV is described. The measurement 
of the cross section is presented together with that of final states containing b quarks. Limits 
on anomalous couplings among neutral gauge bosons are derived. 

2 Data and Monte Carlo Samples 

The data were collected in 1998 by the L3 detector ^ at a/s = 188.7 GeV, and amount to an 
integrated luminosity of 176 pb~^. This energy will be denoted as 189 GeV in the following. 

The EXCALIBUR |1] Monte Carlo is used to generate events belonging to both the signal 
and the background neutral-current four-fermion processes. Background from fermion-pair 
production is described making use of PYTHIA5.72 (e+e" qq(7)), KORALZ4.02 
(e+e~ /i+/x"(7) and e+e" r+r"(7) ) and BHWIDE [Zj (e+e~ e+e~(7)). Background 
from charged-current four-fermion processes is generated with EXCALIBUR for ez/eqq' and 
with i = e,fi,T and K0RALW1.21 jH] for WW production. Contributions from 
multiperipheral processes are modelled by PHOJET 1.05c jH] (e"'"e~ e"'"e~qq) and DIAG36 
[ID] (e+e- ^e+e-^+r). 

The L3 detector response is simulated using the GEANT 3.15 program ^11|, which takes 
into account the effects of energy loss, multiple scattering and showering in the detector. Time 
dependent detector inefficiencies, as measured in the data taking period, are reproduced in 
these simulations. 

The definition of the Z pair signal is unchanged with respect to the generator level phase- 
space cuts of the 183 GeV analysis Those requirements are summarised as follows: the 
invariant mass of both the generated fermion pairs must be between 70 GeV and 105 GeV. 
This criterion has to be satisfied by at least one of the two possible pairings of four same 
flavour fermions. In the case in which fermion pairs can originate from a charged-current 
process (uudd, cess and I'lUei'^i', with i = e, fi, r) the masses of the fermion pairs which could 
come from W decays are required to be either below 75 GeV or above 85 GeV. Events with 
electrons in the final state are rejected if | cos^el > 0.95, where 9e is the electron polar angle. 

The expected cross sections for the different final states are computed with EXCALIBUR. 
A total cross section of 0.662 pb is expected. In this calculation = 0.119 [T^J is included 
for the QCD vertex corrections. The cross section for states with at least one b-quark pair 
amounts to 0.178 pb. 
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3 Event Selection 



All the visible final states of the Z pair decay are investigated, with criteria similar to those 
used at 183 GeV [2,. All selections are based on the identification of two fermion pairs each 
with a mass close to the Z boson mass. The selections are modified to take into account the 
different background composition at the higher ^/s and the changed signal topology due to the 
larger boost of the Z bosons. This boost leads to acollinear and acoplanar fermion pairs. 

3.1 qqi+£- Channel 

A dedicated selection is performed for each of the final states qqe+e~, qqiJ-^fJ.^ and qqr'^r" 
after the application of a common preselection. This requires at least five charged tracks, 15 
calorimetric clusters and a visible energy of more than OA^/s together with two same flavour 
identified leptons. 

Electrons are identified from energy depositions in the electromagnetic calorimeter whose 
shower shape is compatible with those initiated by an electron or a photon. At least one electron 
should have a matched track. Muons are reconstructed from tracks in the muon spectrometer 
pointing to the interaction vertex. Energy depositions in the calorimeters consistent with a 
minimum ionising particle (MIP) which have an associated track are also accepted as second 
muon candidates. For the qqr^r^ channel both a particle-based and a jet-based selection 
are performed. In the first, tau leptons are identified via their decay into isolated electrons or 
muons, or as an isolated low-multiplicity jet with one or three tracks and unit charge. In the 
jet-based selection, the event is forced into four-jets using the DURHAM algorithm. Two 
of the jets must each have less than four tracks. These jets are considered as tau candidates, but 
at least one of them must coincide with a tau candidate defined in the particle-based selection. 

In the electron and the muon channels both the lepton and the jet pair must have an 
opening angle of at least 120°, tightened to 130° for the taus. The invariant mass of the jet- 
jet and the lepton-lepton system after performing a kinematic fit, which imposes energy and 
momentum conservation, must be within 70 GeV and 120 GeV. The events are then subject 
to the DURHAM algorithm requiring lnY34 to be greater than —6.0 for the electron and tau 
channels and —6.5 for the muon one. Y34 is the jet resolution parameter for which the event is 
changed from a four-jet to a three-jet topology. Furthermore, the visible energy in the electron 
channel must be at least 0.8\/s and between 0.6y/s and 0.9^/s for the jet-based tau selection. 

Additional requirements are applied in the tau selection to reduce the radiative qq(7) back- 
ground rejecting events containing a photon of energy larger than 30 GeV. Semileptonic WW 
events are rejected by requiring the transverse missing momentum to be lower than 40 GeV in 
events with no identified electron or muon with energy larger than 40 GeV. 

The kinematic fit is repeated on events that pass at least one of the four selections described 
above with the extra constraint of equal invariant masses for the jet-jet and lepton-lepton 
systems. The distribution of the invariant mass arising from the fit, is shown in Figure 1. 
Table 1 summarises the yield of this selection. 

3.2 qqi/i> Channel 

High multiplicity hadronic events with more than three charged tracks and at least 15 calori- 
metric energy clusters are selected. The event invariant mass must exceed 50 GeV. These 
cuts reduce contributions from purely leptonic two-fermion final states, as well as two-photon 
interactions, while keeping a significant fraction of hadronic events from qq(7) and W-pair 
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Table 1: Data, signal and background Monte Carlo events selected by each analysis and their 
efficiency. The qqi/i/ entries are reported for a selection requirement of 0.5 on the neural network 
output. The t^i'vv figures refer only to electrons and muons. The Monte Carlo statistical 
uncertainties are given on signal and background expectations. 



production. These latter contributions are further reduced by requiring the visible mass to be 
less than 130 GeV and the mass recoiling against the hadronic system to exceed 50 GeV. 

In addition, the transverse momentum is required to be greater than 5 GeV and the lon- 
gitudinal momentum to be smaller than 40% of the visible energy. The energy deposition in 
the forward calorimeters must not exceed 10 GeV and the missing momentum vector must be 
at least 16° away from the beam axis. No electrons, muons or photons with energies above 
20 GeV are allowed in the event and the energy in a 25° azimuthal sector around the missing 
energy direction, £"25, is required to be smaller than 30 GeV. 

A total of 299 events satisfy the selection criteria with 23 and 266 events expected from 
the signal and background Monte Carlo simulations respectively. The dominating background 
is due to chargcd-currcnt four-fermion processes. To differentiate further the ZZ signal from 
the remaining background a neural network is constructed. The inputs to the neural network 
include event shape variables which help to distinguish two-jet from three-jet topologies, the 
sum of invariant and missing masses, the masses of the two jets, the total missing momentum 
and £"25- The use of the neural network increases the signal fraction in the selected sample to 
approximately 60% for large neural network output values, as demonstrated in Figure 2a. The 
efficiency and the yield of this selection are reported in Table 1 for a cut at 0.5, while the full 
spectrum is used for the cross section determination. 

3.3 e.+£-vP Channel 

The selection for £^£^uu is optimised for electron and muon pairs identified as in Section 3.1 
and characterised by an invariant mass, Mu, between 85 GeV and 95 GeV. The requirement 
on the associated track for electrons is dropped and MIPs are not considered. 

In the electron channel only events with a visible energy between 75 GeV and 98 GeV 
are selected; this requirement is loosened to the range 65 GeV — 140 GeV for muons. The 
opening angle of the two electrons must be below 166° and from 143° to 172° for the two 
muons. In order to reduce the background from radiative Bhabha scattering and purely leptonic 
decays of W pairs, the recoil mass to the electron pair is required to be less than 95 GeV. 
The background from other resonant and non-resonant four-fermion processes is reduced by 
performing a kinematic fit imposing the Z mass to the visible pair of Icptons and recalculating 
their four-momenta. The recoil mass, Mrec, after the fit is required to be less than 98 GeV 
for electrons and in excess of 84 GeV but not larger than 98 GeV for muons. The transverse 
momentum has to lie in the range from 4 GeV to 29 GeV. 
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The spectrum of the sum of Mu and Mrec, without applying the kinematic fit, peaked 
around twice the Z mass is presented in Figure 2b. Table 1 lists the efficiencies and the yield 
of the selection. No contribution in the T"^r~ z/t-z/t- signal channel is expected. 

3.4 £+£-£'+£'- Channel 

This selection is based on events with at least four loosely identified leptons of a minimum 
energy of 3 GeV and the subsequent study of just one pair of them. 

First a low multiplicity event preselection is applied, requiring at least two tracks but less 
than 15 calorimetric clusters, with a total visible energy between 0.2^/s and 1.3y^. Elec- 
trons and muons are identified as described in Section 3.1. Low angle electromagnetic showers 
(I cos^l > 0.95) without a matching track are also accepted as electrons. Tau candidates are 
identified as low multiplicity hadronic jets with either one or three tracks in a cone of 10° half 
opening angle. To reject hadronic jets, the energy between 10° and 30° around the tau direction 
must not exceed half of the energy in a cone of 10° half opening angle. To increase the selection 
efficiency, MIPs are also accepted. 

If there are more than four lepton candidates, the four most consistent with energy and 
momentum conservation are chosen. Events are then required to have at least one electron, 
muon, tau or muon-MIP pair. Low angle electrons are not considered in this procedure. If 
more than one such pair is possible, the one with the invariant mass, M«, closest to the Z mass 
is chosen. 

Both Mi£ and the recoil mass, Mj-ec, to this selected lepton pair are required to lie between 
70 GeV and 105 GeV. The data and Monte Carlo distributions for M« + M^ec are shown in 
Figure 2c. Table 1 summarises the total yield of the selection. 

3.5 qqq'q' Channel 

The four-jet selection has to cope with the large QCD and W pair production backgrounds. 
High multiplicity hadronic events are selected by requiring a visible energy between 0.6^^ 
and 1.4-y/i together with parallel and perpendicular imbalances below 0.3^^. Events with an 
identified electron, muon or photon with energy in excess of 65 GeV are discarded. 

The events forced to four jets with the DURHAM algorithm are subjected to a constrained 
fit which rescales the jets to balance momentum while imposing energy conservation. This fit 
reduces greatly the dependence on the calorimeter energy scale. A first neural network [14^ is 
then applied to distinguish events with four genuine quark jets from those with two quark jets 
and two jets from gluon radiation. A cut on the output of this neural network rejects QCD 
background selecting a hadronic sample enhanced in W and Z pairs. 

A second network uses five variables to distinguish Z pairs from W pairs by means of their 
mass difference after the dijet pairings are chosen to minimise the dijet mass difference. The 
variables are the reconstructed dijet mass, the maximum and minimum energy in any jet, the 
average number of charged tracks per jet and the dijet mass difference. A large portion of 
four-jet ZZ decays contains at least one b quark pair, which provides significant distinguishing 
power from W pair decays. A sixth variable, a b-tag discriminant jT^j, is added to the network. 
The output of this network is shown in Figure 2d. A cut on the background enhanced region 
below 0.2 is applied and the rise above 0.8 is due to the b-tagged events. The performances of 
the analysis are summarised in Table 1. 

A simpler sequential analysis is also performed for this channel, taking advantage of the 
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different boost of Z and W pairs and hence investigating just the two dijet opening angles, the 
dijet mass difference and the dijet mean mass. Results compatible with the previous approach 
are obtained limited by a lower purity due to the absence of b-tagging. 

4 Measurement of the ZZ Cross Section 

A binned maximum likelihood fit to each of the variables displayed in Figures 1 and 2 allows 
the determination of the ZZ cross section in the individual final states, as listed in Table 2. 
These are in agreement with the Standard Model values reported in the same table. Assuming 
these predictions as the relative weights of different channels, the ZZ cross section, azz, is found 
to be: 

azz = 0.74l°:i^pb, 

in good agreement with the expected cross section within the signal definition cuts of 0.662 pb. 
The uncertainties are only statistical. 

In the calculation of the cross section the effect of the cross talk between the separate 
channels is found to be negligible. 
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Table 2: Result of the individual cross section fits. 



5 Study of Systematic Errors 

The systematic uncertainties are grouped in correlated and uncorrelated sources among the 
channels. The correlated sources of systematic errors are the background cross sections, the 
LEP energy and the energy scale of the detector. As they modify the shapes of the investigated 
distributions, their effect is evaluated performing a new fit to calculate azz once their values 
are scaled to the extremes listed in Table 3. An uncertainty of 2% is attributed to the measured 
cross section to take into account the difference of the assumed relative weights of the different 
channels with respect to those obtained with the GRC4F ^HI Monte Carlo generator and to 
parametrise the other uncertainties related to their calculation. 

Four sources of systematic uncertainty are uncorrelated among the channels and modify the 
shapes of some of the discriminating distributions. The jet resolution and the charged track 
multiplicity for the qqq'q' selection are scaled as in Table 3. The b-tag procedure and the Monte 
Carlo description of b-hadron jets are taken into account by reweighting the value of the b-tag 
discriminant. Finally the Bhabha Monte Carlo is not sufficient to estimate the background 
to the channel. This estimate is obtained from a fit to the shapes of some selection 

variables whose uncertainty contributes to the total systematic uncertainty. Two additional 
sources of systematic uncertainty are propagated to the total cross section: the Monte Carlo 
statistics and lepton identification. Their uncertainties listed in Table 4 do not affect the shape 
of the discriminating distributions. 

The individual and combined systematic errors are listed in Table 3. 
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The measured cross section is then: 



azz = 0.74+°;}2 (stat.) ±0.04 (syst.) pb. 
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Table 3: Systematic uncertainties on azz and crzz^hhx- 



Channel 


Systematic Source 


Background uncertainty 


Signal uncertainty 


qqi+i' 
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4.7% 


1.7% 
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2.7% 


4.1% 
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Table 4: Sources of uncorrelated systematic uncertainties. 



In terms of the NC02 cross section in which only the two conversion diagrams are considered 
for the double-resonant Z pair production, the cross section reads 

aNC02 _ o.73t°;}l (stat.) ±0.04 (syst.) pb, 

to be compared with a Standard Model expectation of 0.65 pb calculated with the YFSZZ fTfl 
package. 
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6 b Quark Content in ZZ Events 



It is of particular interest to investigate the rate of ZZ events with b quark content. The 
production of the minimal or a supersymmetric Higgs boson would manifest via an enhancement 
of these events and their study complements the dedicated search for such processes The 
expected Standard Model cross section for ZZ^ bbX final states is 0.178 pb. 

The investigation of the ZZ— bbX events proceeds by complementing the analyses of the 
qquu and qqi~^i~ final states described above with a further variable describing the b quark 
content in the event ^3]. Three variables are then considered for each final state: M^c for the 
qqi~^i~ and the neural network output for the qquu analysis together with the b-tag evaluated 
for each of the two hadronic jets. The combination of each of the sets of these three variables 
into a single discriminant proceeds as follows. First the variables are mapped to achieve uniform 
distributions for the background. Then the product of their observed values is calculated event 
by event. Finally the confidence level is calculated for the product of three uniformly distributed 
quantities to be less than the observed product. This confidence level is expected to be low 
for signal and fiat for background. The final discriminant is the negative logarithm of this 
confidence level and is shown in Figure 3. Figure 2d shows the bbqq response from the neural 
network used to select the qqq'q' final states. 

The cross section calculation is performed as above. Results for the individual channels are 
listed in Table 5. The combined result for o"2z_»bbx is: 

^zz^bbx = 0.18l°:°? (stat.) ± 0.02 (syst.) pb. 

This result agrees with the Standard Model expectation and differs from zero at 99.9% confi- 
dence level. In the fit the other ZZ final states are fixed to their Standard Model expectations. 
The systematic uncertainties are evaluated in the same way as for the total cross section and 
are presented in Table 3. Figure 4 displays the measured total and bbX cross sections and their 
expected evolution with a/s. 
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Table 5: Result of the individual ZZ^ bbX cross section fits. 



7 Anomalous Couplings 

A parametrisation of the ZZZ and ZZ7 anomalous couplings is given in Reference [T]. Assuming 
on-shell production of a pair of Z bosons, only four couplings (i = 4,5;V = 7, Z), where 
the V superscript corresponds to an anomalous coupling ZZV, may be different from zero. At 
tree level these couplings are zero in the Standard Model. They are independent from the hf 
couplings that parametrise the possible anomalous ZZ7 vertex P, probed by the e+e~ — > Z7 
process [18 . 

In order to calculate the impact of anomalous couplings on the measured distributions in 
the process e"'"e~ —>■ fff'f', the EXCALIBUR generator is extended ^H]- The matrix elements of 
the Standard Model are supplemented by an additional term containing anomalous couplings, 
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M.acUp'^}, ^, f7) where {p'^} represents the phase space variables and A the hehcities of 
initial and final state fermions. Four-fermion Monte Carlo distributions for non-zero anomalous 
couplings are then obtained by reweighting each event with the factor 



IE (-^4f (K), A) + A^Ac(K}, A, f7)) 



A 



iEl-^4f(K},A) 



X 



where M^fHp'^}, X) is the Standard Model amplitude for the four-fermion final states. An 
average over initial state and a sum over final state helicities are carried out. Initial state 
radiation is taken into account by evaluating the event weight at the centre-of-mass of the 
four-fermion system. Figure 1 displays the effects of an anomalous value of f] obtained by 
reweighting with this technique the four-fermion Monte Carlo events selected by the qqi~^i~ 
analysis. 

The effects of anomalous couplings not only change the ZZ cross section but also the shape 
of the distributions. Using the distributions given in Figures 1 and 2, a binned maximum 
likelihood fit is therefore performed for each of the anomalous couplings f^, fixing the others 
to zero. 

The results of all these fits are compatible with the Standard Model and 95% confidence 
level limits on the couplings are set as follows: 



These limits are still valid for off-shell ZZ production where additional couplings are possible. 
The small asymmetries in these limits are due to the interference term between the anomalous 
coupling diagram and the Standard Model diagrams. Systematic uncertainties on signal and 
background cross sections are taken into account in the derivation of the limits. 

These limits improve by nearly a factor two previously published results at ^/s = 183 GeV 
which are included in the present analysis. 
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Figure 1: Invariant mass after a kinematic fit, Msc, of the lepton pair for the qq£~^i selected 
events. The effect of an anomalous ZZ7 vertex is also shown for a value of its coupling fj = 1.5. 



13 



10 



L3 

ZZ^qqvv 

• Data 

□ ZZMC 

□ Back. MC 



a) 




> 

ev 
c 

^0.5 



L3 

ZZ^IIvv 

• Data 

□ ZZMC 

□ Back. MC 



175 



180 
M„ + M 



rec 



185 
(GeV) 



190 



> 

CJ 
in 
w 
c 

0) 

>0.5 



L3 

ZZ^IIII 

• Data 

□ ZZ MC # # 

□ Back. MC 




" 140 



C) 



M, 



160 
+ M 



180 

rec (GeV) 



200 




Out 



Figure 2: a) Neural network output for qquu selected events, b) Sum of the visible and recoil 
masses for the events selected by the £~^i~uu selection, c) Sum of the invariant and recoil 
masses of the lepton pair closest to Mz for the £'^£~i'^i'~ selected events, d) Output of the 
second neural network for the qqq'q' selection; signal expectations for events with no or at least 
one b quark pair are presented separately. 
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Figure 3: Discriminant variables for a) the hhi~^i and b) the bbz^z^ selections. The last bin 
shows the overflows. 
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Figure 4: Standard Model prediction for the ZZ and ZZ — > bbX cross sections and the corre- 
sponding measurements where statistical and systematic uncertainties are added in quadrature. 
Signal definition cuts implemented with the EXCALIBUR Monte Carlo are applied and a 2% 
uncertainty is associated to the predictions. 
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